% Cholesky decomposition of matrix
% Ref: Egbert & Booker, 1986
% latest date 13.3.2003
   % incomplete

N = 5;
B = rand([1 N*(N+1)/2]);
for I =1:N,
   II=floor((I-1)*I/2);
   fprintf('II = %d\n',II);
            for J=1:I-1,
               JJ = floor((J-1)*J/2);
                  fprintf('JJ+II = %d\n',JJ+II);
           T=B(II+J);
           I1=II;
 
              for K =1:J-1,
              I1=I1+1;
              JJ=JJ+1;
              %            T=T-A(I1)*CONJG(A(JJ))
              % incomplete
           end;
        end;
     end;
     
